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Abstract. In 1991 Tratnik derived two systems of multivariable orthogonal Racah polynomials 
and considered their limit cases. (/-Extensions of these systems are derived, yielding systems of 
multivariable orthogonal (/-Racah polynomials, from which systems of multivariable orthogonal 
(/-Hahn, dual (/-Hahn, (/-Krawtchouk, (/-Meixner, and (/-Charlier polynomials follow as special or 
limit cases. 

Key words and phrases. Multivariable discrete orthogonal polynomials, multivariable basic 
hypergeometric orthogonal polynomials, several variables, multivariable (/-Racah, (/-Hahn, dual 
(/-Hahn, (/-Krawtchouk, (/-Meixner, and (/-Charlier polynomials 

2000 Mathematics Subject Classification. Primary — 33D50; Secondary — 33C50. 
1. Introduction 

The Racah polynomials [1, 2, 25], defined by 
r n (x) = r n (x;a,f3,7,N) 

= (a + l) n (/3 + 7 + l) n (-iV) n 

-n, n + a + (3 + 1, —x, x + 7 — N 



X 4F3 



a + l,/9 + 7 + l,-JV ' X 



(1.1) 



for n = 0, 1, . . . , N, satisfy the discrete orthogonality relation 

N 

J ^2r n (x)r m (x)p(x) = \ n 5 ny m (1-2) 

x=0 

for n,m = 0,1, ... ,N, with the weight function 



p(x) = p(x;a,f3,-y,N) 

1 - N + 2x {"l-N) x {a + l) x {{3 + 1 + l) x {-N) 3 



7 _AT x!( 7 - a- N) x {- 0-N) x ( 7 + l) 3 
and the normalization constant 



;i.3) 



A n = \ n (a,/3,j,N) 

(a + /3 + 2) 7V (-7)jv 



(P + l) N (a-j + l) N 
x (a + + l)n\ (a + l) n {f3 + 1)„(-JV) 
X (a + (3 + l + 2n){a + + l) n 



a - 7 + l) n (P + 7 + l)n (a + P + 2 + N) T 



(1.4) 
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where TV is a nonnegative integer. In 1991 Tratnik [24] extended the Racah polynomials to a system 
of multivariable orthogonal polynomials (in a slightly different notation) 

i?„(x) = # n (x;ai, . . .,a s+1 ,r),N) 

S 

= II r n k ( x k ~ Nfc-i;2iVfc_i + 7? + ak - ai,a fc+ i - l,7V fc _i + a k + x k+1 ,x k+1 - N k _x), (1.5) 
fc=i 

where 

x = (xi, . . . ,x s ), x s+1 =N, n=(n 1 ,...,n s ),N o = 0, 

k k 
N k = y~] n j , 1 < A; < s, afc = ^Oj, 1 < A; < s + 1, (1.6) 
i=i j=i 
and iV s < TV. Clearly, i? n (x) is a polynomial of total degree iV s in the variables y 1 , . . . ,y s with 

y k = x k (x k + a k ), k = l,...,s. 

Tratnik showed that these polynomials satisfy the discrete orthogonality relation 

i?„(x)i? m (x)p(x) = X n 5 n , m (1.7) 

X 

for N S ,M S < N, where M k = Y^j=\ m i an d the summation is over all x = {x\,. . . ,x s ) with 

s 

x k = 0, 1, . . . , N for k = 1, . . . , s, £ njm = fl $n k ,m k -> the weight function is 

k=l 

p(x) = p(x; ai , . . . , a s+1 , rj, N) 

N\T(a s + N + l) { ai ) Xl ( V + l) Xl 



r(a s+ i + N)T(a s+1 + N) Xl l (d - r/k 
A- r(a fc+ i + x fc+ i - x fc )r(a fc+ i + x fc+ i + x fc ) a fc + 2x fc 
^ (xfc+i - x k )\T(a k +x k+1 +x k + 1) a k 



k=l 



and 



A ra («i, • • • , a s+ i,rj, N) 

(a s + TV) jv 3 (r? + 1 - oi - A0jv s 

r(ai - r?)r(a s + /V + l)r(a, +1 - ai + r? + N s + TV + 1) 
X r(a 1 )r(r ? +l)r(a s+1 +iV)r(a 1 -r / + iV) 

s 

II "fc lnk - ^ fc+1 ~ 01 + 71 + Nk + N k~i)n k 
k=l 

T(a k+1 + n k )T(a k - ai + iV fc + 7V fc _i + 1) " 



x 



fl 9) 

T(a fc+1 - oi + v + 2iV fc + 1) J ' V ' 1 

Note that p(x) = if x^+i < for some k < s. He also pointed out the special case of the 
multivariable Hahn polynomials of Karlin and McGregor [14], the limit cases of the multivariable 
Krawtchouk, Meixner and Charlier polynomials, and used permutations of the parameters and 
variables to derive a second system of polynomials that are orthogonal with respect to the weight 
function in (1.8). 



In this paper we derive (/-extensions of Tratnik's i? n (x; ai, . . . , a s+ i, 77, TV) polynomials and of 
his second system of multivariable orthogonal Racah polynomials, from which systems of multi- 
variable orthogonal (/-Halm, dual (/-Halm, (/-Krawtchouk, (/-Meixner, and (/-Charlier polynomials 
follow as special or limit cases. Some (/-extensions of Tratnik's [23] multivariable orthogonal Wilson 
polynomials and of his [21] multivariable biorthogonal generalization of the Wilson polynomials are 
given in our papers [11] and [10], respectively. 

The authors wish to thank Michael Schlosser for his useful comments on the original version 
of this paper. 

2. Multivariable orthogonal (/-Racah polynomials. 

Analogous to (1.1) we define the (/-Racah polynomials [2] by 

r n (x;q) = r n (x;a,b,c,N;q) 

= (a q ,bcq,q- N ;q)n(q N /c) n / 2 

for n = 0, 1, . . . , N, where 



q n ,abq n+1 ,q x ,cq x N 

aq,bcq,q- N ' ,Q,Q 



(2.1) 



n— 1 k 

{a;q)n = I| (1 - aq k ), (a 1: . . .,a k ;q) n = JJ(aj-;g) n , 

fc=0 3=1 

and we use the notation of our book [9]. The power {q N /c) n / 2 is chosen in (2.1) so that certain sym- 
metry properties of the (/-Racah polynomials are satisfied; for example, by the Sears transformation 
formula [9, (2.10.4)] it follows from (2.1) that 

r n (x;a,b,c,N;q) = r n (N - x; b, a, c" 1 , iV; q). (2.2) 

Note that r n (x; q) is a polynomials of degree n in the variable z = q~ x + cq x ~ N . 

Askey and Wilson [2] showed that the (/-Racah polynomials satisfy the orthogonality relation 
[9, (7.2.18)] 

N 

^2 r n (x; q)r m (x; q)p(x; q) = X n (q)S n ,m (2.3) 

x=0 

for n, m = 0, 1, . . . , N, with 

p(x;q) = p(x;a,b,c,N;q) 

1 - cq 2x ~ N (cq~ N ,aq,bcq,q~ N ;q) x 
l-cq- N (q,ca-iq- N ,b-iq- N ,cq;q) x 

and 

A n (?) = X n {a,b,c,N;q) 



{abq)~ x (2.4) 
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(c 1 ,abq 2 ;q) N _ N 

- — — r — (q,aq,bq,aq/c,bcq,q ; 

{aq/c, bq;q) N 

1-abq {abq N+2 ;q) n 



(2.5) 



1 — abq 2n+1 (abq; q) n 
Analogous to (1.5) we define the multivariable g-Racah polynomials by 

R n (-x;q) = i? n (x;ai, . . . , a s+1 ,b, N; q) 

s 

= ]Jr nk (x k - Nk^MkQ 2 ^- 1 /a 1 ,a k+1 q-\A k q x ^ +Nk -\x k+1 - N k _ V ,q), (2.6) 

k=l 

where, in addition to the definitions of x, n, x s+ i and N k given in (1.6), we let 

k 

A = l, A k = Y[ aj , fc = l,...,s + l. (2.7) 

i=i 

Note that -R n (x; q) is a polynomial of total degree N s in the variables z\, . . . , z s where 

z k =q~ Xk +A k q Xk , k = l,...,s. 

In order to derive an orthogonality relation for these polynomials we start by considering the 
following g-extensions of the weight function p(x) in (1.8): 

p(x; q) = p(x; ai , . . . , a s+1 , b, N; q) 

= (q,qA s ;q)N (ai,bq;q) Xl 
(a s+1 ,A a+1 ;q) N (q,ai/b;q) Xl 

„ TT ( <lfc + 1 ' Q)x k+1 -x k (A k +i; q) Xk+1 + Xk {l — A k q 2xk ) 

^ {q;q)x k+1 -x k {qA k ;q) Xk+1+Xk {i - A k ) 

where ci )S , . . . , c S;S are to be determined (see (2.15)) so that the orthogonality relation 

N X s X3 X2 

J2 £ ••• £ £ i?n(x;g)i? m (x;g)p(x;g) = A n (g)5 n , m (2.9) 

a; s =0a; s _i=0 rE2=0a;i=0 

holds for iV s , M s < A and certain normalization constants X a (q), which will be given in (2.14). 

The summation over x\ in (2.9) can be evaluated via (2.3), which is derived in [9, §7.2] by 
means of the terminating 60s summation [9, (2,4,2)], to obtain for s > 2 that 

y> {ai,bq;q) Xl (a 2 ;q)x 2 -x 1 ( A 2;q)x 2 +x 1 { 1 ~ aiq 2xi ) , ,- Xl 

— ' (q,ai/b;q) xi (<i)x 2 -x 1 ( a iq;q)x 2 +x 1 ( 1 - 



X X =0 



x r ni (xx ;b,a 2 q 1 ,a 1 q X2 ,x 2 ;q)r mi (x 1 ;b,a 2 q 1 , cnq X2 , x 2 ; q) 

l-a 2 b {bq,a 2 ,q,A 2 ,bA 2 q/a 1 ;q) ni sbyi n 2 
1 - a 2 bq 2ni (a 2 b;q) ni Vai/ 

* ("0) 

provided we choose ci )S = bq for s > 2. In the single variable case s = 1 it is clear from (2.3) and 
(2.4) that ci,! = 6(7, and hence we conclude that ci jS = 6(7 for s > 1. 



J n k ,m k 



By doing the summations over x\, x 2 , ■ ■ ■ , Xj in (2.9) for small j one is led to conjecture that 
after summing over x\,...,Xj in (2.9) and setting c k>s = a k for 2 < k < s — 1 we have 

yr (l-bA k+1 / ai ) (q,a k+1 ;q) nk (bA k q/ ai ;q) Nk+Nk _ 
11 (1 - bA k+1 q^/ ai ) (bAw/anq^+N^ 

^ (A J+1 ,gb^ +1 /a 1 ;g) J y 3+:ej+1 / b yV-^ v , W'- l// 2 V ' ' \ ' ' (2 n) 

(q,ai/b;q) Xj+1 - N . Vai/ V ai / 

as the sum for j = 1, 2, . . . , s — 1. 

To prove (2.11) by induction on j, assume that j < s — 2 for s > 3, multiply (2.11) by the 
remaining Xj+i-dependent part of the weight function and polynomials, and then sum over Xj+i 
to get 



yr (1 - bA k+1 / ai ) (q,a k+1 ;q) nk (qbA k /a 1 ;q) Nk+Nk _ 
11 (1 - bA k+1 q^/ ai ) (bA^/a^q)^^ 

x (A^ ^ 2 v 9)^+^+1 / Mjg 2jVj+1 \ 

x (oj+2;g)x J . +2 -x J . +1 (A J - +2 ;g)x J . +2 +x J . +1 (l ~ A j+1 q 2x i+i) ^ y Xj+1 
(Q^q)x j+2 -x j+1 (qAj +1 ;q) Xj+2+Xj+1 (l - A j+1 ) 

x r nj+1 (x j+ i - Nf,bA j+1 q 2N > /ai, a j+zq- 1 , A j+1 q N > +x i+ 2 ,x j+2 - Nj-q) 

x r mj+1 (x j+1 - Nj-bAj+iq™ 3 /aiia J+2 q-\A J+1 q N > +x i+\x 3+2 - Nf,q). (2.12) 

Note that the summand above vanishes for < Xj+i < Nj. Hence, the sum is effectively from 
Xj + \ = Nj to Xj+2- Replacing Xj + \ — Nj by x, say, the sum in (2.12) then becomes 

(A j+1 ,qbA j+1 /a 1 ;q) 2 N j («j+2; q) Xj+2 -N j (A j+2 ; q) Xj+2+ N 3 (1 - A J+l q 2N i) 
(q;q)x J+2 -N 3 (qA j+1 ;q) Xj+2+N3 (l - A j+1 ) 

( bA j+1 2Nj+1 y^ Xj+ ^ Ni (1 - A J+1 q 2N ^ 2x ) (A J+1 q 2N >,A ]+2 q^ x ^;q) x 
V ox Q J ^ (l-A j+ig ™>) (q,q N ^-^/a j+2 ;q) x 

v (fe^+x^+Vai.^-^+^^/b^+a 2JV , 



x ^ a l 

x r nj+1 (x;bA j+1 q 2N i /auaj^q' 1 , A j+1 q N ^ +x ^+ 2 , x j+2 - Nj-q) 

x r mj+1 (x;bA j+1 q 2N >/a 1 ,a j+2 q- 1 ,A j+1 q N i +x i+ 2 ,x j+2 - Nj-q) 
_ {qA j +i,qbA j+1 /a 1 ;q) 2Nj (a i+2 ; q) Xj+2 - Nj (A j+2 ;q) Xj+2+Nj 



-bAj+i 2Nj+1 \~ N . 



(<?; q)x j+2 -N j {qAj+i;q) Xj+2+ N j 

) (q, bA j+1 q 2N i +1 / ai ,a j+2 , A j+2 q N ^+ 2 ,q N ^+ 2 ;q) n ^ , 



(1 - bA ]+2 q 2N ^/a 1 ){bA 3+2 q N ^ x ^+ l / ai -q) 



(1 - 6^ +2 ^+^ + v«i)(^+2^7«i;'z)n, +1 
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{q -N 3 -^ /Aj + lMj+2q 2N J+ l /ai . q) ^ +2 _ N 

(bq^-^+ya 1 ,a j+2 ;q) Xj+a . Nj ° n ^> m ^ 

by (2.3). Substituting this into (2.12) and simplifying, we find that the right-hand side of (2.12) is 
(2.11) with j replaced by j + 1. This proves (2.11). 

Now set j = s — 1 in (2.11) and substitute into the left-hand side of (2.9) to find that if we set 
c s s = a s for s > 2, then 



Y R "( x ' 9)^m(x; <?)p(x; q) 

X 

N X s X3 X2 

= Y Y ■ Y Y ^n(x;^)i? m (x ;(7 )p(x;g) 



x s =0 x s _i=0 X2=0 xi=0 
s-1 



(1 - bA k+1 /ai)(q, a k+1 ;q) nk (qbA k /a 1 ;q ) Nk+Nk 

\ k+ i/a 1 ;q) Nk 

{q,A s ,qbA s /a 1 ;q) 2N ^ 1 { A s+iq N ,Q~ N ;q)N s ^ 1 , A „n s ,\-n. 



; =i (1 - M fe+1 ^ /ai)(Mfe+l/ai;g)jVfc+JVfc _ i ™ 

JV _-iV. 



{A s q^\ q ^/a s+1 -q)N^ ^ " ' 



"-Jk- 1 (l-A s q^-^) (A s q^, b -^q™^+\A s+1 q N + N ^,q N ^- N ;q) 3 
^ (1 - A ag ™.-i) (q, ai /6, g i-^.-Vo H U g ^.-+i; g) x 

fr^+i "% Bj {x . bAsq 2N s ^ /aiias+iq -^ As q N + N ^ , TV - 7V s _ i; q) 



x 



ai 

x r ms (x; bArf™'- 1 /^, a a+1 q-\ A s q N+N ^ , iV - iV a _ l5 g) 



A n (g)<5n,m (2.13) 



with 



K(q) = A n (ai, . . . , a s+ i, b, N; q) 

^ {qA s ,qbA s+1 /g 1 ;q) N / ax \^ 
(a s+ i,ai/6;g)jv \g6A s / 

1 -AT / -AT \ 

(q, ak+i 1 q)n k {qbA k /g l ;q) Nk+Nk _ 1 {l - bA k ^j a x ) 
(bA k+1 /a 1 ;q) Nk+Nk _ 1 (l - bA k+1 q 2N >< /ai) 



x (ftA.+ic^+Vai, A.+ig", bq 1 -" / ai ,q- N ; q)i 

n (q,ak+i;q)n k (qbA k /a 1 ;q) Nk+Nk l (l - b^ fc+1/ u, i; 



fc = l 

where, in the weight function ^(x; q) in (2.8) we set 

c\,s = bq for s > 1, Cfc iS = a k for 2 < < s, (2-15) 
which completes the proof of (2.9) with 

(a s+1 ,A s+1 ;q) N (q, ai/6; g) xi Ug/ 
x tt ( a fc+i;g)a; fc+ i-a ;t (^fc+i;g) a;fc+1 +a ;fc (l - ^fcg 23 ''') ^-^ , 2 > 

^ (g;9)x fc+1 -x fc (^fc;g)x fc+1 +x fc (i -A fc ) flfc 

6 



Analogous to Tratnik's [24] permutations of the parameters and variables, we consider the 
following permutations 

ai < — > (A.^)" 1 , a k +i < — ►a s _fc+i, k = 1, . . . , s - 1, 

a s+ i < — > 6g, x fc < — ► AT - a; s _ fc+1 , k = l,...,s. ( 2 -17) 

Clearly the summation region in (2.9) remains unchanged under these permutations, and the weight 
function p(x;q) in (2.16) changes to a multiple of itself, namely, 
(a s+1 ,a 1 /b, A s+1 ;q) N (qbA s \ N 

- f 7 -: x p(x.;au...,a a+1 ,b,N;q). 2.18) 

{a u bq,qA s \q) N \ a x J 

On the other hand the polynomial i? n (x; q) transforms to 

R n ((N — x s ,N — x s -!, . . . ,N — xi); (A s q 2N )~ 1 ,a s ,a s _ 1 , ...,a 2 , bq, a s+ ig _1 , N; q), (2.19) 
which, since (2.6) is not invariant under (2.17), gives a second family of multivariable orthogonal 

s 

(/-Racah polynomials. By replacing n k by n s _ k+ i for k = 1, . . . , s, setting N k = fl nj and using 

j = k 

(2.9) we obtain that the polynomials 
_R n (x; q) = R n (x; a 1 ,..., a s+1 ,b, N; q) 

= r ni (N - N; - x 1 ;A 8+1 ^- 1 /a 1 ,b, q N *~ N /a^N - Nfrq) 

s 

x II r nk (N - N* k+1 - x kj] A 8+1 <? N w- 1 /A k ,a k q-\q N w- N - !C >>-i/A k ,N - N* k+1 - x k - l]q ) 

k=2 

(2.20) 

are a g-extension of Tratnik's [24, (2.12)] second multivariable Racah polynomials and that they 
satisfy the orthogonality relation 

R n (x; q)R m (x; q)p{*; q) = K{q)S n , m (2.21) 

X 

for N s , M s < N, where 

Anfe) = A n (ai, . . . , a s+1 ,b, N; q) 
_ {qAs^qbAs+x/ai.q)^^ / ai \ N 

(a 1 /b,a s+1 ;q) N \qbA s ) 
x (bA s+1 q N+1 / ai , bq 1 ~ N /a 1 , A s+1 q N , q~ N ; q) N . 
(q,qb;q) ni (A s+1 /a 1 ;q) N ; +N *(l - h^+i/ai) 

(bA s+1 /a 1 ;q) N * +N *(l - bA s+l q 2N i / a x ) 
yy (g, Qfc+i; g)n fc+1 (A s+1 /A k+1 ; g)jv; +1 +jy» +2 (1 - A s+1 /gA k ) 
X M (A fl+1 / g A fe ; g ) Jv . +i+JVfc . +2 (l-^ fl+ig 2 ^ + i/ g A fe ) 
3. Some limit cases of (2.9) and (2.21) 

The multivariable g-Racah polynomials in (2.6) and (2.20) contain as limit cases Tratnik's [24] sys- 
tems of multivariable Racah, Hahn, dual Hahn, Krawtchouk, Meixner, and Charlier polynomials. 
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Here we will derive limit cases of the orthogonality relations (2.9) and (2.21) containing multivari- 
able g-Hahn, dual g-Hahn, g-Krawtchouk, g-Meixner, and g-Charlier polynomials. We start with 
the g-analogue of the dual Hahn polynomials defined by 

d n (x; b, c, N; q) = lim r n (x; a, b, c, N; q) 



/ \ n / 2 

= (bcq,q~ N ;q)n[q N /cj 3^2 
for n = 0, 1, . . . , N, and let 

D n (x;q) = D n (x;ai, . . . ,a s+1 ,N;q) 

= lim R n (x.;ai,...,a a+ i,b,N;< 



-n g x^x N 

bcq,q- N 



(3.1) 



s 



Y[d nk (x k - jVjfc_i;o fc+ i9 1 ,A k q Xk + 1+Nk - 1 ,x k+1 -iV fc _i;<?) 



(3.2) 



fc=i 



where N s < N and x s+ \ = N. Then, by taking the 6^0 limit of (2.9) we get the orthogonality 
relation 

^ Ai(x; q)D m (x; g)pn(x; q) = A D (n; q)5 ri , m (3.3) 

for N s , M s < N, were 

p D (x;q) = pu(x;ai, . . .,a a+1 ,N;q) 
(q,qA s ;q) N (a x ;g) 
(a s+ i, A s+ i;g)iv 



n 

fc=i 



(Qfc + i;g)a: fc + i-a: fc (^fc + i;g)x fc+ i+a; fc (l ~ A k q Xk ) 

(qA k ;q) Xk+1+Xk (l - A k ) 



l k > 



(3.4) 



A D (n;g) = A D (n;ai, . . . , a a+ i,N; q) 

(qA s ; q) N (A s+1 q N , q~ N ; q) Na 



(-qA s )- N q-(2) 



X 



(«s+i; q)n 

s 

Y[(q,a k+1 ;q) nk , (3.5) 
fc=i 

and the summation in (3.3) is over the same region as in (2.13). 

On the other hand, as b — > the orthogonality relation (2.21) approaches the limit 

^ T n (x; q)T m (x; q)p T (*; q) = A T (n; q)5 n , m (3.6) 

X 

for N s , M s < N, where 

T n (x;g) = T n (x;ai, . . .,a 8+1 ,N;q) 

= d ni (N - N; - x 1 ;A s+1 q N+N ^-\q N ^ N /a 1 ,N- N^;q) 

1 

-,7V- N k+1 - Xk-nq), 

(3.7) 



II rnu(N ~ N* k+1 - x k ; ^q^~\ * S^l 

A k q 



N'-N-Xk-! 



k=2 



Ai 



pr(x; q) = />d(x; q) is the same weight function as in (3.4), and 
A T (n;g) = A T (n;ai,. . .,a s+1 ,N;q) 

-(-^ S )-V (2) 



( a s+i; q)n 

x (^s+i?^, 9 _iV ; ?)iv* {q; q) ni (A s+1 /ai;q) N * +N * 



s 1 (g,a fc+ i;g) njb+1 (A s+1 /A k+1 ;q) N . +i+N . +a (l - A s+1 /qA k ) 
x 11 



fc=i 

If we multiply (2.9) by 



-"fc /Mi. A/f _i_ 



ai 



n 

fc=i 1 

and take the limit b — > 00, we obtain the limiting relation 

X] D n( x ; 9)^m( x ; <?W ( x , 9) = ^D* (n; q)5 xlyI 

X 

for 7V S , M s < N, where 

D n(*;q) = D^(x; ai ,...,a s+1 ,N;q) 



- [I <C (** " N k _ 1 ;a k+1 q- 1 ,A k q x ^ +N ^ , x k+l - N k _ l]q ) 



k=i 



with 



d* n (x;b,c,N;q) = lira (aq) n r n (x;a,b,c,N;q) 

a— >oo 

= (bcq,q- N ;qU-irq( n 2)(q N /c)^ 2 3 
for n = 0,1,... ,iV, 

PD*(x;g) = pd*(x;oi, . . . ,a s+ i,iV;g) 



Tb ry X S~*f1"^ 

1 y ! °y 



-A' 



rc+1 



(g^A,;^ (ai;< 



(a s+ i, A s+ i;g)jv 

(A k+1 ;q) Xk+1+Xk (l-A k q 2 ^) 



n 

fc=i 



(ajfe+i;g) 



(9;9)x fc+1 -x fc (?-4fc;g)a; fc + 1 +a; fc (l ~ -4* 



and 



\ D *(n;q) = A D *(n;ai, . . .,a s+1 ,N;q) 

= ^^(-a s+i rq^A^^ 
{a s +i;q)N 



x l[(q,a k+1 ;q) nk (A k q N ^ 1 )- 2nk Q 



N k _ 1 + l\-2n k „N k _ 1 -N k -N k -N k _ 1 



"k+1 



k=l 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



(3.13) 



Note that pr>*(-x.;q) = {a\/q) Xl ) q^^> po(x;g) and that D*(x;g) is not a constant (independent 
of x) multiple of L> n (x;g). Similarly, by inverting the base q in the orthogonality relation (3.9) it 
follows that (3.9) is equivalent to (3.3). 



To obtain g-analogues of the multivariable Hahn polynomials in Karlin and McGregor [14] and 
in Tratnik [24], multiply (2.9) by a^ B , take the limit a\ — > 0, replace b by a\ and a k by qa k for 
k = 2, 3, . . . , s + 1, and make the change of variables y\ = x\, y k = x k — Xk-i for k = 2, 3, . . . , s. 
This yields the multivariable (7-Hahn polynomial orthogonality relation 

5]ffn(y;9)ffm(y;?K(y;?) = AH(n ;(? )5 n>m (3.14) 

y 

for iV s , M s < AT, where 

H n (x;q) = iJ n (x;ai, . . . , a s+ i, TV; g) 

= n^( y *-^- i ; A *« 2JVfc+fc_1 ' a *+ i ' r *+ i - iV *- i '«) ( 3 - 15 ) 

fc=l 

with the single-variable g-Hahn polynomials defined by 

h n (x; a, b, N; q) = lim(cq~ N ) n/2 r n (x; a, b, c, N; q) 



= (aq,Q N ;q)n 302 

for n = 0,1,... ,N, 



q n ,abq n+1 ,q x 
aq,q 



(3.16) 



pH(y;q) = pHiy;^,. . . ,a s+1 ,N;q) 

(q~ N ;q)Y s - Ys rr (i a ^q) 



- (q-N /as+1 . q)Y3 a ^Ll. {q . q)yk 
\ H (n;q) = A#(n;ai, . . .,a s+1 ,N;q) 

= (A s+1 q s+1 ;q) N +N s (q~ N ;q)N s ( ^ gJV „ +s)) - jV( ^N 3q ( N 2 °) 

(qa s+1 ;q) N 

TT (g» q a k+i,q)n k (A k g k ; q) Nk+Nk _ 1 (1 - A k+l q k ) / fc+2iVfc 

M OW^W+^a - ,W+™*) l^ 9 J ' (3 - 18) 

Yfc = ^ yj for 1 < /c < s so that Y\= y\= x\ and x k = Y k tor I < k < s, Y s+ i = x s+ i = N, and 

.7 = 1 

the summation is over all y with j/^ = 0,1,..., and Y" s < AT". Also see Dunkl [4] and Rosengren [18]. 

Clearly the summation region in (3.14) and the weight function in (3.17) are invariant under any 
permutation of the labels (1, 2, . . . , s). If we set y s +i = N — Y s , then they are also invariant (apart 
from the renormalization of the weight function) under any permutation of the labels (l,2,...,s+l); 
i.e., under any simultaneous permutation of (ai, a^, ■ ■ ■ , a s +i) an d (yi, yz, ■ ■ ■ , y s +i)- Since the 
polynomials in (3.15) are generally not invariant under these permutations, this generates distinct 
systems satisfying the orthogonality relation (3.14) via these permutations, which are ^-extensions 
of those mentioned at the top of page 2341 in Tratnik [24] . 
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JJ up^ ((zafc) -v fcj (317) 



Next, consider the g-Krawtchouk polynomials denned by 

k n (x; b, N; q) = lim (aq)~ n h n (x; a, b, N; q) 

a—too 

= (q~ N ; q)n(-l) n q^ 2^1 (q~ n , g"*! g"^! q, bq n+1 ) (3.19) 

for n = 0, 1, . . . , N. In view of (3.19), if we multiply (3.14) by a~[ 2Na , let a\ — > 00, and then replace 
Ofc by afc_i for 2 < < s + 1, we obtain the multivariable g-Krawtchouk orthogonality relation 

2 K n (y; g)# m (y; gW(y; g) = A^(n; g)<5 n>m (3.20) 
y 

for N s , M s < N, where 

K n (y,q) =K n (y;ai,...,a a ,N;q) 

S 

= Y[k nj (Y j -N j _ 1 ;a j ,Y j+1 -N j _ 1 ;q), (3.21) 
3=1 



pK(y,q) =p K (y;au...,a a ,N;q) 



(q- N ;q) 



Y, 



.(-l)vi q -( V i)aJ Y > 



(q;q)yAq N /a s ;q) Ya 

xn^,-i)- y ', (3-22) 



Ax(n; g) = A^(n; ai, . . . , a s , N; q) 

= to-""}*' (A s _ ig ^)-^(^^ 1 )^(-I)^ g (^)+( JV+ 2 JVS ) 

xn^^;^^,!^-)-^^^" 1 ?-^ (3-23) 

and the summation is over the same region as in (3.14). This orthogonality relation can also be 
derived by starting with the a\ — > limit case of (3.14). 

Now consider the g-Meixner polynomials defined by 

M n {q-*; a, c; q) = 2 fa (q~ n , q~ x ; aq; q, -q n+1 /c) (3.24) 

for n = 0, 1, . . ., which satisfy the orthogonality relation 



£ M n (q- X ; a, c; g)M m (g"*; a, c; q) <? q ® 
— ' (q, -acq; q) x 

q- n 6 n , m - (3.25) 



x=0 (q,-acq;q) 
(-c;g)oo(g, -g/c;g) 



(-acg;g)oo(ag;g) n 
See [9, Ex. 7.12] and [15, §3.13]. 
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To find multivariable orthogonal g-Meixner polynomials we, formally, replace the upper limit 
of summation N for the sum over x s in (2.9) by oo and replace q~ N by an arbitrary number /3, 
say. Then, from (2.16), the factor 

(q, gA s ;q) N (a s+1 ;q) N - Xs (A s+1 ;q) N+Xs 
(a s+1 ,A s+1 ;q) N (q; q) N - Xe {qA s ; q) N+Xe 
(q- N ,A s+1 q N ;q) Xa 



(qi- N /a s+1 ,A s q^;q) Xs 

gets replaced by 



(q/a s+1 ) x % 



and, from (2.14), the factor 



(p,A a+1 /0;g) Xs 
((3q/a s+1 ,qA s /(3;q) Xa s_t 

(qA s ,qbA s+1 /g 1 ;q) N / ai ■ • 
(a s+ i,ai/b;q) N \qbA s 



{qA s ,a x q N /bA s+1 ;q) 



N 



gets replaced by 



(q 1 N /a s+1 ,a 1 /b;q) N 
(qA s ,a 1 (3/bA s+1 ,q/a s+1 ,a 1 /b/3; q) c 



((3q/a s+1 ,a 1 /b, qA s /(3, ai/bA a+1 ; q)^ 

via the identity (a;g)jv = (a; q)oo/ iflq N \ q)oo in [9, (1-5)]. Hence, by multiplying both sides of 
(2.9) by {bq/a\^ 2 )~ Ns ~ Ms , taking the limits b — > oo and ai — > 0, and replacing at by afc_i for 
2 < fe < s + 1, we are led to the limit case 



OO X s 



Yj Yj ■■■ Y M n {x;q)M m {yL;q)p M (x;q) 

i s =0 i s _i=0 x 1 =0 

= \ M (n;q)S n , m , (3.26) 



where 



X n (x; q) = M n (yi; a u . . . ,a s , j3;q) 

~s-l 

Y[{q Nk -'- Xk+1 ]q) nk {-l) nk q^)+ nkNk ^Al k J? 
Lfc=i 

xM nk {q Nk --- Xk -q Nk --- Xk +-~\ -q/a k ; q)] 
x ((3q N ^-q) ne (-l)^q( n 2) +n ^A^_{ 2 

x M ns (q Ns - 1 - Xs ^q Ns - 1 -\-q/a s ;q), (3.27) 
Px(x;g) = p^(x;ai, . . .,a s ,(3;q) 



(qP/a s ;q) a 



Va s _ia s / (g;g) Xl 
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with ao = 1, and 



X M (n;q) = Ax(n;oi, . . . ,a s , (3;q) 

(q/a s ;q) 



, W;q) Ns q N ^~ 1) A^ 

x fl(q,a k ;q) nk a- k Nk - Nk -\ N ^- N *. (3.29) 

k=i 

This orthogonality relation can be verified by starting with the g-Meixner orthogonality relation 
(3.25) and proceeding as in the derivation of (2.9). 

The (3 — > limit of (3.26) gives the following multivariable extension of the g-Charlier polyno- 
mial orthogonality 

OO X s X2 

^2 $Z • • • X] C n(x;<?)C m (x;<7)/? c (x;g) 

i s =0 i s _i=0 xi=0 

= \c(n;q)8 n , m , (3.30) 

where 

C n (x;<?) = C„(x;ai, . .. ,a s ;q) 

"s-1 

Y[{q N ^- x ^-q) nk {-\) n *qC2)+^ N ^A n k k _l? 
_fc=i 

x M ni (q N K-i-** ; - g / 0fc ; g)] 

x (-l)"s q C 2 e )+n s N s _ 1A n L /2 

xc ne (q N °-i- x °;-q/a s ;q), (3.31) 
Pc(x;g) = p c (x;ai, . . . ,a s ;<?) 

= (^-y t^- q d) ft 1 ^k±^ an , ( 3.32) 

Ac(n;g) = A c (n;ai, ... ,a s ;g) 

= (9/a.;g)oo^ r ' (JV '- 1) ^ fffe a k ; ^a^^- 1 ^^, (3.33) 

fc=i 

with ao = 1 and the o-Charlier polynomial defined by 



c n (q x ;a;q) = lim M n (q x ;f3,a;q) 



= 2Mq- n ,Q- X ;0;q,-q n+i /a). (3.34) 

See [9, Ex. 7.13] and [15, §3.23] for the orthogonality relation and other properties of the g-Charlier 
polynomials. Of course (3.3), (3.6), (3.9), (3.14), (3.20), and (3.30) can also be derived directly 
from their single variable special cases by proceeding as in the derivation of (2.9) and (2.21). 
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Other discrete multivariable extensions of the Racah and Hahn polynomials are considered in 
[3], [5], [13], [17], [18], [20], [22] and [26]. For some related results see, e.g., [6], [7], [8], [9], [12], 
[16] and [19]. In view of the nonnegativity results for the linearization coefficients and for kernels 
containing products of (/-Racah polynomials in [7] and [8] and the resulting convolution structure 
and positive summability methods, it would be of interest to see if any of these nonnegativity results 
can be extended to some of the multivariable orthogonal polynomials considered in this paper. 
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